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Abstract 

The Standard Model flavor structure can be explained in theories where the fermions are localized 
on different points in a compact extra dimension. We show that models with two bulk scalars 
compactified on an orbifold can produce such separations in a natural way. We study the shapes 
and overlaps of the fermion wave functions. We show that, generically, realistic models of Gaussian 
overlaps are unnatural since they require very large Yukawa couplings between the fermions and 
the bulk scalars. We give an example of a five dimensional two scalar model that accounts naturally 
for the observed quark masses, mixing angles and CP violation. 
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I. INTRODUCTION 



The Standard Model (SM) flavor puzzle is to understand the origin of the smallness 
and hierarchy of most of the fermion masses and mixing angles. The likely solution of this 
puzzle requires the existence of a fundamental theory that generates the observed patterns 
in a natural way, namely without small dimensionless parameters. One such framework is 
due to Arkani-Hamed and Schmaltz (AS) 0. The idea is to separate the various SM fermion 
fields inside compact extra dimensions. In that case the four dimensional Yukawa couplings 
between two fermions are suppressed by the overlap of their zero mode wave functions. Once 
the widths of the wave functions are smaller than their separation in the extra dimensions, 
small and hierarchical Yukawa couplings in four dimensions are induced. 

Various scenarios based on this idea had been constructed. In some models, only the 
fermion wave functions have non trivial shapes in the extra dimension EL S @' S B @» 01 • m 
other cases, also the Higss field is used to generate the flavor structure ||, |8|, ||. A specific 
model in five dimensions (5d) that generates the observed fermion masses and mixing angles 
was found by Mirabelli and Schmaltz (MS) ||. This model, however, cannot accommodate 
the observed CP violation in the kaon and B systems ||10| . Experimental tests of the split 
fermion idea were studied in fl|, pTT | . 



While the AS framework is very attractive there is no complete realization of it. In |T]| 
a model with one infinite extra dimension, x 5 , was considered where domain wall fermions 
were used to localize the fermions and to separate them. The basic idea is that there is a 
scalar, $, which develops an x 5 dependent vacuum expectation value (vev), h(x 5 ). Then, 
the fermions have £5 dependent effective masses m e ff = A/i(xs) — M, where M is a bare 
mass term and A is the Yukawa coupling between the fermion and The minimum of 
the right handed (RH) zero mode and the maximum of the left handed (LH) zero modes is 
where the effective mass vanishes, namely at 2™ such that A/i(x™) = M. The LH chirality 
state, with the maximum at is localized around x™. The RH mode, however, is localized 
at infinity, namely it is non-normalizable and it drops out of the spectrum. The separation 
between the LH modes is achieved by introducing a different M for the different fields. 

The generalization of this original realization of the split fermions idea to the realistic case 
of a compact dimension is not trivial. The problem is that in a compact space both LH and 
RH zero modes are normalizable, and thus the model is not really chiral. Another problem, 
which is related to the first one, is that in a compact space the domain wall configuration is 
not stable. 

Some ideas toward solving these problems had been proposed by Georgi, Grant and Hailu 
(GGH) || and by Kaplan and Tait (KT) |4]]. They considered 5d scenarios with one extra 
dimension compactified on an S\jZ 2 orbifold. The Z 2 orbifold symmetry projects out one 
chiral zero mode, making this a chiral theory in four dimensions (4d). A scalar, $, which 
is odd under the Z 2 symmetry is also introduced. Since $ is odd, its vev must vanish on 
the orbifold fixed points. When $ develops a non zero vev in the bulk, an x§ dependent 
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vev is generated. Fermion mass terms of the from are allowed while bare mass 

terms, are forbidden by the Z 2 symmetry. Therefore, the zero modes are localized 

around one of the boundaries since these are the only places where their effective masses 
vanish. The sign of the Yukawa coupling between the fermion and the scalar determines on 
which boundary the fermion is localized. The idea of split fermions where each fermion is 
stuck at a different point in the bulk cannot be achieved. The zero modes can be localized 
only around one of the orbifold fixed points. An additional flavor structure in the Yukawa 
couplings is needed in order to naturally generate the observed fermion parameters in such 
two location models j3|. More complicated models, using large non renormalizable 1 terms 
H or boundary terms M can localize fermions in the bulk. 

The problems of this setup can be traced down to two points. The first is that the scalar 
vev vanishes only on the boundaries. Thus, the fermions cannot be localized within the 
bulk. The second is that mass terms are absent. In the original AS model the scalar was 
used to localize the fermion zero modes, and different constant mass terms were used to split 
them. In the simple orbifold models, while the scalar can be used to localize the fermions, 
there is no simple mechanism that splits them. 

We are therefore interested in improving upon the models of [§,[§]. We propose to use 
two scalar fields that couple to the fermions. This rather minor modification of the GGH 
and KT scenarios can produce localization in the bulk and at the same time does not affect 
the appealing features of these models. The advantage that two scalar models have over 
one scalar models is that the effective mass can vanish in the bulk. Intuitively the picture 
is as follows. For one scalar the sign of the Yukawa coupling determines the boundary 
where the fermion is localized |J. Once a second scalar with opposite sign Yukawa coupling 
is introduced, the picture is more complicated. The second scalar tends to localize the 
fermion on the other boundary. In some cases one scalar is dominant and localization is 
at the boundary determined by the sign of its Yukawa coupling and the second scalar only 
modifies the shape of the wave function. In other cases, however, the tension between the 
two scalar results in a compromise: A configuration where the fermion is localized in the 
bulk. 

In section II we introduce the two scalar model and investigate its scalar sector. In section 
III we study the shapes of the fermion wave functions in various limits of the model and 
their implications for model building. In section IV we present a two scalar model which 
accounts naturally for the quark masses, mixing angles and CP violation. We conclude in 
section V. 



More precisely, we refer here to 5d terms which after dimensional reduction produce non renormalizable 
operators in the 4d effective theory. 
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II. TWO SCALAR MODEL 



We now move to study the features of the scalar sector of our model. We explain how it 
can generate localization and separation of fermion wave functions in the extra dimension. 

The space-time of our model is described by ordinary 4d space-time and an additional 
space dimension compactified on an S1/Z2 orbifold. The physical region is defined as < 
x 5 < L and x 5 = 0, L are the orbifold fixed points. The model includes two scalar fields, $j 
(i = 1,2), and a fermion field, \E'. The scalars and are odd while ^ l is even under the 
Z2 orbifold symmetry. For simplicity, we assume that there is no interaction between the 
two scalars. We further consider the case where the two scalars develop £5 dependent vevs. 
The conditions for this case to be realized can be found by a straightforward generalization 
of the single scalar case (see e.g. 0). Here we only remark that having an x 5 dependent 
vev is a generic situation. Then, the action is 



~2\2 



where fi and Aj are real with Aj > 0. (Our notation is such that capital letters run over the 
five dimensions while Greek letters over the four dimensions.) 

It is convenient to use dimensionless parameters. We therefore use rescaled fields 



j0_ 

VI' 



and introduce the following dimensionless quantities 



v 



(2) 



u 



X5 
L ' 



2\i ViL, 



fi 



fi 



X 



Q2 j 2 



(3) 



Note that < u < 1. We order the two scalars such that a± < a 2 and further define / = f\ 
and a = a±. Then, the action can be written as 



S 



dx„ / du 



1 1 a 2 

^-d^id^ - -d 5 (pid 5 ipi - -j- (ip 2 - \y 



(4) 



where d 5 = 

We move to investigate the functions hi(u) = ((pi){u). For one scalar they were investi- 
gated in detail in @, 0. Since we assume that there is no interaction between the scalars 
we can use the one scalar model results for each of the scalars. Neglecting the Yukawa 
interactions, and working in the limit where [0] 



di > 1, 



(5) 
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hi(u) 



(7) 



we have 

hi(u) = tanh [ajii] tanh [o»(l — «)] • (6) 

Since hi(u) is symmetric under u — » 1 — w we consider only the < w < 1/2 range. We will 
make use of the following approximations 

' 1 for it ^> 1 /a* , 
v a>iU for it <C 1 /dj . 

Of particular interest is the following function 

^(«)=Ai(«)-X/i 2 («). (8) 

Given the values of tij and recalling that a x < a 2 , the sign of g(u) is determined by X as 
follows: 

(i) X > 1 



sign [0 («)] = -1 



(9) 



(ii) X < a 1 /a 2 



sign[flf(«)] = 1 



(10) 



(hi) 1 > X > ai/a2 



sign 



u< — 
a 2 



sign 



(11) 



Eq. (|Tl|) teaches us that when condition (iii) is satisfied g(u) vanishes not only on the 
boundaries. Namely, there exists a point < u max < 1/2 such that g(u max ) = 0. (Note that 
since hi(u) is symmetric under u — > 1 — u if g{u) crosses zero in the bulk, it crosses it twice, 
at itmax and 1 — it max .) 

While we cannot solve analytically for w max in the general case, we can do it in what we 
denote as the constant mass approximation 



ai/a 2 1 



In this approximation w max is given by 



arctanh(X) X ^,^q n 

Mmax W ^ = — + 0(X 3 ) 

a a 



(12) 



(13) 



Figs, (lb) and (lc) show the function g{u) (the thin black curves) in the cases described 
by eq. (|TTD for typical values of parameters. The figures demonstrate the main point of this 
section: The two scalar model yields an effective scalar profile which has several zeros. In 
particular, one of them is not an orbifold fixed point. Thus, we can construct models in 
which the fermions are localized at various points in the fifth dimension and not only at the 
orbifold fixed points. 
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III. THE FERMION ZERO MODE 



In the following we analyze the functional behavior of the fermion zero mode wave function 
and, once more fermions are introduced, the overlaps between their zero mode wave functions 
in the fifth dimension. Though we do it in the context of a two scalar model we shall see 
that some of our conclusions are rather generic. They apply to other realizations of the split 
fermions scenario. 

Following the standard treatment (see e.g. [|1|]) we find the differential equation for y(u), 
the fifth dimension wave function of ^> l : 

1 dy(u) . . „. 

, x V ' =-fag(u), 14 
y(u) ou 

where g(u) is defined in (pi). The solution is given by 



y(u) = N exp 



—fa / g(w) dw 



(15) 



where N is the normalization factor. The local maxima of y(u) are at points in which g{u) 
vanishes. Which of these maxima is the global maximum depends on the model parameters, 
and can be at any of these points. 



A. A single scalar model 

We begin our analysis by setting X to zero in eq. (§D and returning to the single scalar 
case |l], ||. Then, the fermion zero mode is given by 

y(u) ~ iV exp {— / In [cosh(cra)]} for -u<l/2. (16) 

We shall analyze y(u) in two limits, large and small /. We show that these are related to 
the two scenarios studied by AS [l|] and KT || respectively. We assume that / is positive 
and thus the maximum of y(u) is at u = 0. Then, we have 



y(u) fa N exp 



fa 2 



u 2 



for u<l/o; (17) 



2 

y(u) w N exp [-f(au - In 2)] for 1/a <u< 1/2 . (18) 

When eq. (|T7|) is applicable, the fermion wave function is approximately a Gaussian with a 
width 

1 



(19) 



When eq. fll8 ) is applicable, the fermion wave function is approximately an exponential 
with a width 

r~-L. (20) 
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Substituting u = 1/a in fllTD we learn that when / 3> 1 the exponential part is negligible 
and then is of a Gaussian shape. Similarly, substituting u — 1/a in flTBp we learn that 
for / < 1 the Gaussian part is negligible. 



B. A two scalar model 

We now move to the two scalar scenario. While the integral in flT5|) can be solved an- 
alytically, the result is not very illuminating. We therefore consider y(u) in various limits 
related to the three cases of eqs. (|D|), flTU| ) and (|TT|). 

When X > 1 or X < g^/g^ the function g(u) vanishes only on the fixed points and y(u) 
is localized at either u — or u — 1. This case is similar to the single scalar model (with 
a zero mass) discussed above. The corresponding fermion wave function will be of the form 
of either Gaussian or exponent as shown in eqs. ( |T7D or (|18|). In the limit of large and 
with / of order unity it reproduces the KT exponential model. 

When 1 > X > ai/ci2 the fermions can be localized in the bulk. It is useful to discuss this 
case in the constant mass approximation [see eq. (|T2"D]. In this approximation the second 
term in g(u) can be treated constant and then 

y(u) ~ N exp {-/ [In [cosh(aw)] - X}} for u < 1/2 . (21) 



Substituting X = ML/af in (21) , with M being the 5d fermion mass, the model reproduces 
the constant mass models of [jl], 0. Thus all our following analysis and conclusions apply to 
these models as well. 

We like to approximate y(u) near its maxima. The situation in this case is more involved 
since y(u) has two maxima and which is the dominant one depends on the interplay between 
the values of X, a and /. Inspecting Eq. ( |I5|) we see that this is determined by the relative 
size of I p and I n defined as 



1-ttr, 



I p = g(u)du, I n = g{u)du. (22) 

JO J Umax 

Once I p > I n the global maximum is at one of the fixed points, while in the opposite case 
it is at w m ax or 1 — M max - We define Xf, to be the transition point, namely, the point where 
Ip = In- Using the constant mass approximation [eq. (|21|)1 and the linear approximation for 
hi [eq. (0)] we find the following estimate for the transition point: 

Xb ~ ^—^ . (23) 
b 2a -3 V ; 

We therefore consider three cases: 

• X > Xfe. The global maximum is at one of the orbifold fixed points. This case is 
similar to the case of a function with only one maximum and does not lead to a new 
structure. 
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• X < Xt,. The dominant maximum is in the bulk, < M max < 1, and y(u) is approxi- 
mately given by 

y(u) ~ N exp 

This is a Gaussian with a maximum at w max and a width 



ay/f{l - X*) " 

For X <C 1 the width is independent of X [as in eq. fll9l) 1, and the maximum location 
is linear in X. This is the case discussed in |], [2[] (with the substitution X = ML/af). 
The generic situation, however, is that the peak of the wave function is not linear with 
X and that its width grows with the distance of the peak from the origin. 

• X ~ Xf,. The two maxima are significant. This case cannot be described analytically 
but it is a combination of the above two cases. We remark that this transition region 
between the other two cases considered above becomes sharper as the value of / 
increases. Since we are working with moderate values of / it does not require fine 
tuning of X to have a situation where both peaks are important. 

In figures (la) -(lc) we plot the effective two scalar vev, g(u), and the corresponding 
fermion wave function for these three cases. 



fa 2 (1-X 2 



[u - u r 



for 



\U — U r 



< 1/a 



(24) 



C. Wave functions overlaps 

We extend the model by including the SM fermion fields and one Higgs field. After 
compactification, the zero modes of these fields consist of the 4d SM fields. In particular, we 
are interested in how the 4d Yukawa interactions are generated form the 5d theory. Consider, 
for example, the 5d Yukawa interaction between a quark doublet (Q), a down type quark 
singlet (D) and the Higgs field (H) 

£ = -J±=QHD + h.c., (26) 

where M* is the natural scale of the 5d theory. After compactification the 4d Yukawa 
interaction is given by 

C = Y 5 Kd L hd R + h.c, (27) 
where lower case letters represent 4d zero mode fields. The overlap K is defined as 

K= [ y D (u)y Q (u)du, (28) 
Jo 

where yo (vq) is the fifth dimension wave function of the singlet (doublet) field. We assumed 
here that the Higgs vev is flat in u. While this is not generally the case, we assume that its 
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shape is not far from flat, and thus approximating it by a flat profile is reasonable. Assuming 
y 5 ~ 0(1) the flavor structure arises from the various overlap integrals JTJ. 

Suppose that we consider a realistic 5d model which accounts for the quark masses and 
mixing angles by Gaussian like suppressions. (Such overlaps were assumed, for example, 
in constructing the MS model and in many other papers see e.g. 0, ||.) The model is 
specified by the Yukawa couplings /* and effective masses X % of the nine SM quark fields. 
We ask what are the ranges of the 5d flavor parameters in order for the above model to yield 
the correct masses and CKM matrix elements. 

In order to obtain Gaussian suppressed overlaps and to reproduce the MS model several 
conditions are required. First, the significant overlap between any two fermion wave function 
has to be in the region where both wave functions are, to a good approximation, Gaussian. 
Second, their maxima have to scale linearly with X. The first condition is satisfied in the 
large / limit. The second condition requires small X. 

To make our discussion concrete we assume that the most distant fermions have X < 0.3 
which implies "Umax ~ 0.3/a. Using Eq. ( p4f) we find that the overlap between the Gaussian 
wave functions, assuming universal couplings /, is given by exp(— fX 2 ). In || it was found 
that a viable model of quark masses and mixing requires a separation of roughly 18 standard 
deviations between the two most distant quarks. Then, for X < 0.3 we find 

/>2(-^Y~ 10 4 . (29) 



0.3 

Assuming perturbativity, this implies that A, [introduced in (|])] must be much smaller than 
their natural values JT2|. Note that to a large extent, condition ( p9[ ) is independent of the 
fundamental theory parameters such as L, a and M#. Thus, it is generic that models with 
Gaussian overlaps require very large couplings /. 

While the model of pure Gaussian wave functions seems unnatural, using the analysis 
of subsection |111 B| it is possible to construct models which account naturally for the quark 



flavor parameters. This is demonstrated in the next section. 



IV. TOWARD A REALISTIC NATURAL MODEL 

As already mentioned, several models that produce the observed fermion masses and 
mixing angles had been constructed. Some of them use Higgs wave functions that are not 
flat to generate the flavor structure. We do not discuss these models here. Instead we discuss 
models that assume, similar to our model, that the Higgs vev do not play a significant role 
in generating the observed flavor patterns. In the Appendix we discuss the GGH model 
since our model is not directly related to it. Below we describe the KT model M and our 
improved version of it. 

The KT setup is as follows. LH fermions are localized on one of the orbifold fixed points, 
and RH fermions on the other one. The different fermion wave functions have different 
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widths such that the lighter the fermion is, its LH and RH wave functions are narrower. 
This induces hierarchical structures for the corresponding mass matrices. Yet, in order to 
produce the observed quark masses and mixing angles, this structure is not enough and the 
5d Yukawa couplings must also have some flavor structure 0] . 

One reason for that, for example, is related to the large top mass. In order to generate 
the 0(1) 4d Yukawa coupling, the overlap between y® 3 and y U3 has to be large. Since these 
wave functions are localized on different orbifold fixed points, it implies that at least one of 
them has to be very broad. Both possibilities, however are problematic. When y® 3 is broad, 
a large mass for the b quark is generated, in contrast to observation. When y m is broad the 
overlaps between it and the other doublets are large. Then, large CKM mixing angles are 
generated, in disagreement with the data. Similar problems are encountered when one tries 
to account for the flavor parameters of the other generations. 

The KT model is significantly improved once some of the fields are localized in the bulk. 
For example, the problem with the third generation is solved once y® 3 and y U3 are localized 
near each other far from any of the orbifold fixed points. Then, their overlap is large and 
large m t is produced. At the same time, the overlaps of y U3 with the other quark doublet 
wave functions are small, thus producing small CKM mixing angles. The b quark mass is 
produced by constructing y d3 to have a double peak shape, such that it has moderate overlap 
with y® 3 . 

Based on these ideas we construct an example of a configuration which roughly reproduces 
the correct quark masses and mixing angles. We choose the following parameters: 



J I 



d\ = 4 . 



"1,2,3 
1 



-1,1,28 



-1,-1,21 



a 2 = 12 ; 

X Ql > 2 ' 3 = -4,-1.9,0.73; 
x «i,2,3 = -2.3,2.3,0.78; 
X dl ' 2 ' 3 = -2,-1.4,0.87. 



(30) 



The resultant shapes of the various wave functions are shown in figs. (2a) and (2b) for the 
up and down type quarks respectively. The mass matrices of the up type and down type 
quarks are given, up to order one coefficients, by (in MeV units) 



/ 1.9 3.1 • 10 2 
4.7 7.2 • 10 2 
\ 1.7 - 10 2 3.8 -10 3 



6.4 • 10 2 ' 
4.2 • 10 3 
1.6 • 10 5 



M d ~ 



/ 5.3 4.0 -10 1 
1.3 • 10 1 9.3 -10 1 
V2.4-10 2 7.2 ■ 10 2 



5.3 • 10 1 

3.4 • 10 3 



(31) 



Their eigenvalues are given by 



m u ~ 1.9 MeV. 
m d ~ 5.3 MeV . 



m c ~ 7.2 • 10 2 MeV 
m, ~ 9.3 • 10 1 MeV . 



m t ~ 1.6 ■ 10 5 MeV, 
m b ~ 3.4 • 10 3 MeV . 



(32) 



which agree with the known quark masses, calculated at the Mz scale |L3| , |14[| . The rotation 
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matrices from the left that diagonalize the mass matrices are given by 

1 2A 0.7 • A 3 \ / 1 2A 0.4 -A 3 ' 



VI ~ I 2A 1 0.5 • A 2 ] Vt - 
,0.4A 3 0.6 -A 2 1 



2A 1 0.4 -A 2 | , (33) 

V0.5-A 3 0.4 -A 2 1 



where A = 0.22. The absolute values of the elements of the CKM matrix, Vqkm = {Yl)^li 
are given to leading order in powers of A by 



1 


A 


A 3 \ 


A 


1 


A 2 


A 3 


A 2 


1 ) 



(34) 



The values of the quark mixing angles agree with the experimental data [TJJ. Assuming 
arbitrary phases we find for the Jarlskog measure of CP violation pj 



J~10~ 5 , (35) 



in agreement with the value obtained from global fit to the available data |T3| . We conclude 
that the above model also reproduces the observed size of CP violation. (The same is true 
also in the original KT model |4]].) This is in contrast to the MS model, where the observed 
order of magnitude of J cannot be accommodated [|H]] . The reason is that in the MS model 
there are very small entries due to the strong Gaussian suppressions. Here this is not the 
case. The matrices M u and M d do not contain such tiny entries. Actually, the generic 
situation is that split fermion models in five dimensions that produce the correct masses 
and mixing angles also account for the observed CP violation. 

The example above is only a demonstration that the observed fermion masses and mixing 
angles can be generated. The apparent small deviations from the experimental values are 
due to the unknown Yukawa couplings. The point we emphasize is that the flavor hierarchy 
arises only from the overlap between the different wave functions. The unknown 0(1) 
Yukawa couplings do not carry any hierarchical flavor structure. 



V. DISCUSSION AND CONCLUSION 



There are several related issues that we did not study and seem to us worth further 
investigation. We did not discuss the lepton sector. It is easy to generate a model for the 
charged lepton masses. The issue of neutrino masses is more complicated and may require 
additional ingredients. We also did not investigate the Higgs vev shape. While we assumed 
that it is flat, generally this is not the case. It seems to us that it can be flat enough in order 
not to upset our general conclusions. However, it is interesting to fully understand the vev 
shape and its implications. 

Our model can be extended, for example, by adding more scalars. Then, the number 
of the fermion wave functions maxima can be more then two, up to the number of scalars. 
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One of the maxima is an orbifold fixed point, and the others are in the bulk. Having more 
than one maxima in the bulk can be used to obtain a richer structure. However, this extra 
flexibility is not needed in order to generate the observed quark masses and mixing angles. 

The framework we considered cannot be directly related to the fundamental theory of 
nature. It must be a low energy limit of a more fundamental theory that should also explain 
the compactification and stabilization of the orbifold setup. Therefore, in general, non 
renormalizable terms and terms localized on the orbifold fixed points can be present in the 
effective theory. Since such terms modify the shape of the fermion wave functions they were 
used in || and Q to generate the needed structure. These terms, however, have to be very 
large in order to have significant effects. In our case, however, we do not need them. Thus, 
assuming that they are at their natural size, they are not expected to significantly affect the 
important ingredients of our model: The number of maxima, their widths and their rough 
locations. 

So far we considered only the zero modes. Heavy modes, namely, the Kaluza-Klein tower 
of the SM fields and heavy SM singlets fields, cannot be ignored since they can significantly 
contribute to SM rare processes. In particular, proton decay, N—N oscillation, meson mixing 
and lepton number violating processes. The most severe bound on the properties of the heavy 
fields is usually obtained form proton decay data. In models where the fundamental scale 
is low, very small overlap between the leptons and quarks wave functions was proposed as 
a way to avoid rapid proton decay QXJ] . The required tiny overlap can be achieved in models 
with Gaussian wave functions. However, as we explained, such models require significant 
fine tuning. Thus, proton stability can be explained by having a very high fundamental 
scale or by adding new ingredients to the model. On the other hand, the data on other rare 
processes can be accommodated in models with a relatively low fundamental scale (see e.g. 

§)■ 

To conclude we summarize our main results. We showed that realistic models of split 
fermions can be constructed using two scalars in a world with one compact extra dimension 
compactified on an orbifold. We found that models of split fermions with Gaussian overlaps 
require large 5d Yukawa couplings, and are therefore unnatural. We construct a realistic 
two scalar model where the overlap is not purely Gaussian which accounts naturally for the 
quark flavor parameters. 
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APPENDIX A: THE GGH MODEL 



Here we analyze the GGH mechanism of fermion separation 0. We show that in this 
model it is rather unnatural to obtain sharp localization and separation of the fermions at 
the same time. 

With our scaling [see eqs. (0) and @], the GGH model is described by 



dx, 



dulip 



fa 



1 



L 2 M 2 



1> 



d 2 



L -d^tpd^i - -d 5 ipd 5 ip - — [y - 1 



(Al) 



When a is large, g(u) is approximated by 
g(u) ~ /atanh(au 



2/x i 9 / \ 

1 — r^-sech iau) 
b 1 



for u<l/2 



(A2) 



with 6 = LM*/a ^> 1. The fermion wave function in the fifth dimension is given by [see eq. 

In [cosh (cm)] H — ^-secli 2 (au) j . (A3) 



(0)] 

y(u< l/2)« iVexpj-/ 
The maximum is found at 

Mmax = -arctanhi/l 
a V 

and the width is given approximately by 

1 1 



fe2 

2ft 



(A4) 



(A5) 



a^f] tanh(aw max ) ' 
Note few problematic points: 

• In order to have a maximum in the bulk, very large coupling is needed, ft > b 2 /2 ||. 

• When ft ^> b 2 /2, the peak of the wave function is far from the origin. In that case, 
however, w m ax oc log(ft). Consequently, to separate the fermions in the extra dimension, 
different flavors must have large hierarchy in their corresponding values of ft. 

• When M max is small the width is proportional to 1/M max - Consequently, a very large / 
is needed in order to get a sharply localized wave functions. 



[1] N. Arkani-Hamed and M. Schmaltz, Phys. Rev. D 61, 033005 (2000) |hep-ph/9903417 |. 



13 



[2] E. A. Mirabelli and M. Schmaltz, Phys. Rev. D 61, 113011 (2000) jhep-ph/9912265H . 

[3] H. Georgi, A. K. Grant and G. Hailu, Phys. Rev. D 63, 064027 (2001) |iep-ph/0007350 |. 

[4] D. E. Kaplan and T. M. Tait, JHEP 0111, 051 (2001) |hep-ph/0110T26C . 

[5] H. V. Klapdor-Kleingrothaus and U. Sarkar, Phys. Lett. B 541, 332 (2002) [|hep-ph/0201226| ; 



P. Q. Hung and M. Seco, |hep-ph/0111013| ; Y. Uehara, JHEP 0112, 034 (2001) [jEep- 



ph/010729^ ]; G. Barenboim, G. C. Branco, A. de Gouvea and M. N. Rebelo, Phys. Rev. 



D 64, 073005 (2001) ||hep-ph/01043l2l ; J. L. Crooks, J. O. Dunn and P. H. Frampton, Astro- 

phys. J. 546, LI (2001) |astro-ph/0002089i . 
[6] T. Matsuda, Phys. Rev. D 66, 047301 (2002) [|hep-ph/020533l] ; Phys. Rev. D 66, 023508 

(2002) ||hep-ph/0204307|| ; Phys. Rev. D 65, 107302 (2002) [|hep-ph/0202258 |. 
[7] G. Perez, [hep-ph/02081~02l . 

[8] N. Haba and N. Maru, |hep-ph/020406g| ; |hep-ph/0202196| ;Phys. Lett. B 532, 93 (2002) [ |hep 



ph/0201216H ; N. Haba, N. Maru and N. Nakamura, |hep-ph/0209009| ; N. Maru, Phys. Lett. 
B 522, 117 (2001) lhep-ph/0108002|1 ; M. Kakizaki and M. Yamaguchi, ftiep-ph/0110~266 
M. Maru, N. Sakai, Y. Sakamura and R. Sugisaka, Nucl. Phys. B 616, 47 (2001) {iep- 



th/010720^1 ; M. Kakizaki and M. Yamaguchi, Prog. Theor. Phys. 107, 433 (2002) [|hcp- 



| ph/0104103H . 

[9] D. E. Kaplan and T. M. Tait, JHEP 0006, 020 (2000) [|hep-ph/0004200[| ; G. R. Dvali and 
M. A. Shifman, Phys. Lett. B 475, 295 (2000) [|hep-ph/000T072l ; A. Hebecker and J. March- 
Russell, Phys. Lett. B 541, 338 (2002) |iep-ph/0205143| ; F. Del Aguila and J. Santiago, JHEP 
0203, 010 (2002) [ fcrXiv:hep-ph/0111047| j. 
[10] G. C. Branco, A. de Gouvea and M. N. Rebelo, Phys. Lett. B 506, 115 (2001) [[hep 



ph/001228GH . 



[11] N. Arkani-Hamed, Y. Grossman and M. Schmaltz, Phys. Rev. D 61, 115004 (2000) [pcp- 



ph/990941lH ; T. Han, G. D. Kribs and B. McElrath, |hep-ph/0207003t W. F. Chang, I. L. Ho 



and J. N. Ng, |hep-ph/0203~2~l1| S. Nussinov and R. Shrock, Phys. Lett. B 526, 137 (2002) 
[ arXiv:hep-ph/0101340|l ; Phys. Rev. Lett. 88, 171601 (2002) [|hep-ph/0112337|] ; D. J. Chung 
and T. Dent, Phys. Rev. D 66, 023501 (2002) [|hep-ph/0112360|l ; G. Barenboim, G. C. Branco, 
A. de Gouvea and M. N. Rebelo, Phys. Rev. D 64, 073005 (2001) Hhep-ph/0104312l ; 
T. G. Rizzo, Phys. Rev. D 64, 015003 (2001) [|hep-ph/0101278| ; A. Masiero, M. Peloso, 
L. Sorbo and R. Tabbash, Phys. Rev. D 62, 063515 (2000) [|hep-ph/0003312|1 ; A. Delgado, 
A. Pomarol and M. Quiros, JHEP 0001, 030 (2000) lhep-ph/9911252|1 . 

[12] Z. Chacko, M. A. Luty and E. Ponton, JHEP 0007, 036 (2000) [>iep-ph/9909248| . 

[13] K. Hagiwara et al. [Particle Data Group Collaboration], Phys. Rev. D 66, 010001 (2002). 

[14] See e.g.: H. Fusaoka and Y. Koide, Phys. Rev. D 57, 3986 (1998) flhep-ph/9712201|] ; M. 
Carena et. ai, hep-ph/0010338| . 

[15] C. Jarlskog, Phys. Rev. Lett. 55, 1039 (1985). 



14 




(la) (lb) (lc) 

Fig. 1: The typical shapes of the scalar vevs and the fermion wave functions. The thin black 
curves correspond to the effective scalars vev, g(u) [eq. @]. The thick gray curves correspond to 
the fermion wave functions, y(u) [eq. (^jj)]- The relevant parameters are a± = 5, 02 = 17, and 
X = 1.15,0.75,0.85 for figs, (la), (lb) and (lc) respectively. In this example X^ ~ 0.87 [eq. 
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(2a) (2b) 

Fig. 2: The configuration of the quark wave functions in the extra dimension. Figs (2a), (2b) 
are related to overlaps between the up and down type quarks respectively, y^ 1 - 2 * 3 are the dotted, 
dashed, and dashed-dotted curves respectively (shown in both figs). y U1 ' 2 ^ and y dl ' 2 > 3 are the thin 
black, mid gray and thick light curves of figs. (2a) and (2b) respectively. The relevant parameters 
are given in [eq. 
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